A magnetic field diagnostic for sonoluminescence 
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This study is motivated by the extraordinary process of single bubble sonoluminescence (SBSL), 
where an acoustically driven spherical shock is thought to power the emitted radiation. We propose 
new experiments using an external magnetic field which can induce anisotropies in both the shock 
propagation and radiation pattern. The effects will depend on the temperature, conductivity, and 
size of the radiating region. Our predictions suggest that such a laboratory experiment could serve 

- - , as an important diagnostic in placing bounds on these parameters and understanding the physics of 
\^ • sonoluminescence. 

Os ■ 
G^ ■ 

Q Sonoluminescence, first discovered in 1933 is a remarkable phenomenon whereby sound is converted into light. 

1^ Recently, the process involving a single small gas bubble trapped in a degassed liquid such as water, which is then 
Q ' acoustically driven, has been studied Since the sound wavelength is much larger than the bubble radius R{t), 
] (see Fig. |l|), the bubble feels a uniform pressure which varies with time. The gas bubble then undergoes complicated 
,—1 nonlinear oscillations. At a threshold of the driving pressure of about 1.15 atm, short (< 50 ps), intense pulses of 
light are radiated. The acoustic energy of the system is thus focused into a small region within the bubble for very 
' short times. Recent interest in sonoluminesence has been in part motivated by possible technological applications in 
1^ [ biophysics, sonochemistry and reactions at ultra-high temperatures 0|. 

■ Experiments have probed the effects of ambient temperature, liquid composition, gas composition, driving intensity, 
I and frequency on the phenomenon. In the experiments of Ref. |3|, there appears to be a sharp transition from a 

■ nonradiating oscillating bubble to a luminescing one. Spectral fits to blackbody or Bremsstrahlung emission are 
inconclusive but suggest temperatures of at least 5000°K 

Theoretical attention has focused on bubble dynamics, usually employing variations of the Rayleigh-Plesset equa- 
tion to model the bubble motion and stability Pfl. Furthermore, studies imply that as the gas bubble shrinks, 
' it launches an inwardly propagating shock wave M. The shock collapses to the center, rebounds, and can hit the 
^ ' liquid/gas bubble wall from which it was launched. Radiation is thought to be emitted immediately after the shock 
I rebounds from the center, where extremely high temperatures and ionization are predicted |^ . Numerical simulations 
■ O suggest temperatures of 10* °K much higher than those implied by experiments. Thermal blackbody and/or 
^ Bremsstrahlung radiation, Casimir effects, or decay of excited molecular states, have all been suggested as relevant 
^ processes involved in the luminescence, but the emission spectra in SBSL are so featureless that distinguishing among 
!• radiation mechanisms is difficult. Despite attention on the effects of bubble dynamics, diffusion, and material compo- 
. ^ [ sition on sonoluminescence, our understanding of the underlying physical mechanisms remains an unsorted mixture 

■ of nonlinear hydrodynamics, shock physics and chemical reaction kinetics. However, one emerging criterion for SBSL 
5h I is the necessity of stability of spherical symmetry in bubble structure R{t) [|| . The transient symmetric stability of 

- - ■ the converging shock is also though to be crucial for SBSL. 

In this Letter, we suggest the use of an externally applied magnetic field as a diagnostic by catalyzing the breaking 
of spherical stability and putting bounds on the temperature, conductivity, and ionization of the sonoluminescing 
bubble. The external magnetic field can break spherical symmetry of the shock, bubble, or radiation patterns and 
disrupt or alter sonoluminescent behavior. 

Effects of B on bubble structure - We assume that the bubble gas is ionized at some point in the acoustic cycle. 
Adiabatic compression of the gas bubble R{t) is probably sufficient to partially ionize the gas; Saha's equation for 
temperatures of 7000°K give roughly 5% ionization. First consider the possibility that during the adiabatic collapse, 
a spherical region of plasma has a high conductivity Up such that the magnetic Reynolds number Rm = (JpvR* 3> 1 
(where are typical velocities and the radius of the ionized region). Alfven's theorem of flux freezing holds, and 

the application of B^xt adds an magnetic stiffness in the direction perpendicular to B^xt- The ionized region then 
collapses asymmetrically, with the poles compressing more than the equator. If Rm 2> 1 is not achieved by initial 
adiabatic compression, shocks launched by the liquid/gas interface may further compress, heat, and ionize the plasma, 
consistent with numerical simulations 

Shocks form when the hydrodynamic velocity exceeds a characteristic group velocity. In the presence of a magnetic 
field, characteristic group velocities depend on the magnitude and orientation of the applied field. For Rm ^ 1, the 
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MHD modes have group velocities 



(1) 



and V = va cosO, where v\ = /(Anp) and Cg is the hydrodynamic sound speed {B is the magnetic field in the region 
of interest and p is the plasma density) . These are the magneto-acoustic waves and the Alfven wave respectively. The 
angular dependence of these velocities will break the spherical symmetry of a plasma shock. 

We consider a shock just after rebound since here the upstream region was previously shocked during pre-rebound 
and more likely to be in a high Rm state. We will first assume Rm ^ 1 both upstream and downstream. (The 
possibly important but rather involved case of ionizing shocks will not be treated here.) The governing equations 
are the mass, momentum, and energy conservation, and the steady state Maxwell's equations and Ohm's law, E ~ 
c/ (47rcr)V X B — (y/c) x B, vahd for = mec/eBo < collision time (c, me, e are the speed of light, electron mass, 
and electron charge respectively) . Locally integrating the conservation laws perpendicular to a shock, these equations 
are V • B = V • (pw) = and, pi 



r]n X V X B = fix {v x B) 
pv {h-Vvi + ^UiV -v) = {B-n)Bi/A'K — pvi{v-n) ~ PMni 
^n-{\/xB)xB ^ \fi-v{u+ ^pv^ + Pm) ~ ^{B-v)(B-h) 
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where [S']g = S{x) — S{Q) and u — P/(7 — l)p is the internal energy of the gas (P, 77, and 7 — Cp/C^ are the 
pressure, shear viscosity and adiabatic index of the gas, respectively). The limits straddling the shock are far outside 
the shock bubble (0), where all quantities are uniform, and x, a normal distance downstream. The total pressure 
Pm = P + B^/Stt (for Rm > 1, the E"^ contribution is smaller than B"^ by (w/c)^, which is small; for Rm ^ 1, E^ is 
smaller by (v / c)"^ / Rm) ■ Equations include dissipation due to plasma shear viscosity v and resistivity rj = / (47rcr). 
We can neglect the contribution of the any short intense radiation burst to the electromagnetic stress tensor if we 
focus only on the shock conditions just before or after the radiative burst. Compression ratios and pressure jumps 
depend on 0, the angle between the Bq and the shock normal. We thus expect a nonuniform MHD shock forming 
first at the poles where 9 = Q and i?o offers no additional stiffness. 

In the high conductivity {Rm — ^ 00) limit, all quantities vary over a thin shock front such that the fields at a; = 0^ 
are uniform and gradient terms vanish when we consider the quantities away from the sharp transition region. Upon 
solving the resulting equations (Rankinc-Hugoniot) wc obtain, 

a [7M2(1 - l/s) - A + 1] = (a2 _ 1) sin^ e 
^ (A - 1) + - 1) = 4/3-^1 - ") sin^ e - 7M2 tan^ (l - f 

where a = {(3jAIq — 2 cos^ d)/{(3^MQ / s — 2 cos^ 9), (3 = SttPq / Bq is the ratio of the particle to magnetic pressure, and 
Mq — uq/cs is the Mach number. The roots of Eqs. (H), s and A (> 1), determine the thermodynamically allowed 
shock states. 

Fig. H shows the local shock structure for 7 = 5/3 and at varying Mq and (3. Except for low {3 and Mq, the 
compression and density ratios, A and s, incur their greatest jumps at the poles. Numerically predicted pressures are 
on the order 10^° dyne/cm^ thus, fields of Bq ~ lOT are required for (3 w 0{1). The angle 9 can more or less be 
associated with the angle of B^xt in Fig. |l| Deformations of the ionized region before shock formation will stretch the 
x-axes in Fig. || depending on f3. Fig. ||, however, shows the difference in jumps between 9 = 0, Tr/2. For large (3, the 
shock will be nearly spherical, but magnetic field compression near the equator would tend to decrease 6 throughout 
the hemisphere. When (3 is small, the shock surface will be oblate owing to the lateral magnetic stiffness, and again, 
smaller 9 will approximate most of the surface. In fact, the shock solution ends at a switch-off shock for low enough 
f3. Dissipative processes will smear this region. 

In the limit (low conductivity, or small velocity or length scales) the magnetic field and the plasma fiow 

are decoupled, and the plasma can slide freely through the field. The discontinuities in s, A and v approach those of 
a hydrodynamic shock and the effect of the external i?-field for Rm ^ 1 is expected to be small. Though the result 
given by (||) is correct for all Rm, for small Rm the i?- field dissipation length L becomes much larger than the viscous 
dissipation length S ~ the scale over which pressure and velocity vary. 

Instead of solving the nonlinear differential equations (|^) we simply estimate the variation of B across S at 
9 = tt/2, where the largest effect occurs. Most of the variation in v, p, and P will occur over S; far enough away 
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from the shock the differences in the uniform states (0 and 1) are given by (^). Using dB/dx ~ (Bq — Bi)/L, and 
evaluating quantities in Eqs. ph across S, s and A to order Rm/P are 



' - ' (M„^(7rl)+7)^ [— ' 



(4) 



where s = (7 + l)Af^/(7 + Mq(7 — 1)) and A — (4A'/q — 7 + l)/(7 + 1) are are the density and compression ratios 
of a pure hydrodynamic shock. The angular dependences obey /(O) = g{0) = and /(±7r/2) = g(±7r/2) = 1. As 
expected, the shock anisotropy is small if Rm //? is small. 

Effects of B on radiation - Even when B^xt does not alter the shock structure, it can have other effects on SBSL. 
We now consider the effects on the SBSL radiation. Shock simulations suggest maximum temperatures of ~ 10^ °K 
while spectral fits from radiation emanating from the bubble center suggest a more modest ~ 10'*°K The 
average thermal velocity at these temperatures is < 0.1c, hence, we will only consider the nonrelativistic conditions 
for observing anisotropic cyclotron radiation in competition with isotropic Bremsstrahlung or blackbody emission. 
We ignore collective effects and thus require ujb > (^p, the plasma frequency. For appropriately dense plasmas, 
LUp « 10" - lO^Vs, requiring Bq > 1 - lOT. 

If shock asymmetry is negligible {e.g. if Rm <C 1), then the plasma feels the externally applied magnetic field Bq ~ 
Bextz and charge trajectories obey Ve = (e/me)ve x Bq, where Ve is the electron velocity. The dipole approximation 
for cyclotron emission power per unit solid angle of a collection of uncorrelated electrons is given by [0 

dL^/dO.'^ {nee'^/4TTc^){n x {n x Ve))^ = (riee'^/47rm^c^)t;iB^(l - sin^ 6'sin2 93) (5) 

where Ue is the ionized electron density, n is in the direction of the observer, 9 is the angle between the field and the 
line of sight, and ip is the azimuthal angle to be averaged over. Note that for <C 50ps, also requires Bq > I — lOT. 

To get the integrated power as a function of 6, we average over the uncorrelated Ve using f{ve,x,t)d^Ve — 
exp[— + v'^)]v±dv±dvzdip, where vt = {"^ksT /meY^'^ ■ We also impose a cutoff on the velocity, Vz,v± < £cWb, 
where £c = min{A, R*}, A = AkgT^ / riee'^ , is the mean free path for electron collisions, and R* is the size of the ionized 
region. This ensures that the charge can coherently twist around the magnetic field sufficiently to radiate near ojb- 
The cyclotron power is then 

L,{9) = n,{e^B^vl/87Tmy) [l - (A + 1)6"^] Erf(\/A)(l + cos^ 9), (6) 

where A = {£cUJB)'^m/2kBT. When 4 = A < i?c, A cx T^, but when 4 = < A, A cx 1/T. Thus the maximum range 
of v± , Vz contributing to would occur where A ~ R^. When A > 1, 1^(9) = n^{e-^v^B^/16Trmlc^){l + cos 6), the 
standard expression |]. When A < 1, 

1,(9) ^ neK^r)-^/' (^Q elcolil + cos^ 9) + 0(A^/2). (7) 

To determine when anisotropic emission may be observable, we compare the cyclotron intensities with those from 
isotropic Bremsstrahlung and blackbody emissions. We consider a singly ionized species in a neutral plasma, justified 
because the Debye length <ti R*. Near ujb, the isotropic Bremsstrahlung radiation power in the dipole approximation 
is. 



where &mm ~ 4e^/7rmet;|i is a minimum impact parameter, and where T is the cyclotron line width or the nar- 
row bandwidth of a detector. Detector bandwidths can be made very narrow (< lcm~^), so we use F ~ ojcoi ^ 
{2k})T /uif.X'^Y/'^ ^ 7 X 10^^(ne/10^^)(r/10^K)^'^/^s~^ from coUisional broadening, which dominates Doppler broad- 
ening for T < lO^K. (We have assumed the neutral plasma to have a density ~ lO^kg/m"^ consisting primarily of singly 
ionized nitrogen.) The power radiated by a blackbody at low frequencies (for B = lOT, ujb is in the IR/microwave 
region ~ 1.8 x lO^Vs) is Lbb{ujB) =i TujlkBT/ (tt^c^) - IO^^ub/I.S x 10"sec-i)2(T/106K)(F/7 x 10i2sec-i)erg/sec. 
Similarly, Lbrem 6.4 x 10^'^{ne/lO^'^cm-^)^{T/10^K)-'^/'^(T/7 x 10i2sgc-i)erg/sec, where the Gaunt factor loga- 
rithm in (ph is ~ 1. Now A ~ 10"'', for these characteristic parameters, and Lc ~ 6 x 10^^(ne/10^^cm~'^)(a;B/1.3 x 



3 



10-^^sec~^)''(4/10~^cm)^(T/10^K)~^/^erg/sec so will dominate. In this estimate we have taken 4 ^ R* ■ Note 
that the dominance of Lc depends strongly on ■ 

Discussion - We have illustrated how a magnetic field i?ext may break the symmetry of a collapsing ionized region 
or a propagating shock (Fig. ^ in two limits of Rm. Using Spitzer's formula to estimate Cp at 10^ °K and typical 
velocity and length scales from shock simulations we conservatively estimate Rm ~ 10~^. However, in general, 
consider Guderly's solution for a spherical hydrodynamic shock: r{t) oc t", where t = is the time the imploding 
shock collapses. Here, v{t) oc t"~^, and Rm — crpt^^~^. For noninteracting gases, n > 1/2, implying Rm — + as the 
shock converges {t — > 0) and rebounds. It is possible that n < 1/2 for compressed materials with large 7 — Cp/Cy. 
Experiments using gases that have rt < 1/2 when compressed and heated, if they exist, would be in the i?M 00 
regime, where shocks will be strongly perturbed by B^xt (Fig. ||). 

In either case, large or small Rm, anisotropics develop during the evolution of a shock. These anisotropics are 
greater for smaller f3 (larger Bq), but vanish as Rm / P 0. Spherically converging hydrodynamic shocks are inherently 
unstable. If luminescence requires a spherically converging shock, Bext may further destabilize the shock leading to 
reduced heating upon implosion and destroy luminescence. Comparison of experiments with numerical results can 
then produce an estimate of Rm. Furthermore, the 0-dcpendent pressure exerted on the liquid/gas interface by a 
recolliding shock may destabilize the nonlinear bubble oscillations. An external magnetic field can thus indirectly 
alter the region of bubble oscillation stability |8|. It is thought that bubble stability is crucial for SBSL. Experiments 
with external fields can determine the relative importance of symmetric shocks for radiation phenomenon and stable 
bubble oscillations in SBSL. Differences between multiple bubble sonoluminescence (MBSL) and SBSL can also be 
probed; if only SBSL depends strongly on symmetric bubble oscillations, an external field would affect SBSL and not 
MBSL. 

For low Rm, especially with strong magnetic fields, the radiation can be expected to have an anisotropic cyclotron 
component given by (|^). Experimental determination of the luminescing spectrum is interrupted by the surrounding 
water which cuts out light with wavelength below ^ 220 nm. We therefore suggest that structure may be observed 
at longer IR and microwave wavelengths, near ujb and its higher harmonics. An angular dependence of the form 
~ 1 + cos^ 6 may be observed if enough signal can be collected from a narrow bandwidth detector set at nws . Given 
the stable, oscillatory nature of SBSL, this seems feasible. One may be able to discern structure at these harmonics by 
comparing the complicated low frequency spectra with and without Bq. It is straightforward to estimate the cyclotron 
contribution by comparing Lc/{Lc + Lhrem + ^bb) at = and 7r/2 using (§) and (|). If no effects are observed with 
strong Bext, the interpretation depends on the optical depth: If the plasma were optically thick, blackbody would be 
the only emission mode. If the plasma were optically thin, the absence of cyclotron emission means that the velocity 
cutoff determined by A would be very small because ujcoi 3> (j^s- Lhrem would therefore dominate. Since Lhrem and 
Lhb differ by many orders of magnitude near ldb, the two cases should be distinguishable. Radiation anisotropy may 
be detected at low Rm- 

The spectroscopy of a sonoluminescing bubble in a magnetic field can help put bounds on relevant physi- 
cal parameters such as charge density, and temperature. In particular, Le oc T~^/^ when ic — R* < A, but 
Lc oc T^^/^ when £c — X < R* ■ The temperature when R* — A, above which Lc oc T'^/"^, is given by 
Tcrit ~ 10^(i?c/10~^cm)^/^(ne/10^^cm~^)^/^K. Since increasing the acoustic driving pressure increases T, the func- 
tional change in Le (anisotropy) could help constrain T . Finally, note that experiments hitherto have time-integrated 
the radiation. A time dependent spectral analysis of asymmetric radiation near ojb can be a means of monitoring 
the temperature and charge density of the heat-up and cool-down phases of the bubble oscillation, in addition to the 
50ps burst. Experiments of this nature may help clarify the mechanisms involved in sonolumincscncc. 

T.C. thanks M. Brenner and R.V.E. Lovelace for helpful discussions and N. Kannan for assistance in plotting Fig. 
2(a). 
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FIG. 1. The sonoluminescing bubble in an external magnetic field Bqz. The liquid/gas interface and the shock front are 
denoted by R{t) and r{t) respectively. The external sound field supplies the pressure Pa- 



FIG. 2. The pressure and density ratios for Rm 3> 1. (a) A as a function of 9 for Mo = 1.5 (lower 3 thin curves), and 
Mo = 3 (upper 3 thick curves). Solid, dotted, and dashed curves correspond to /3 = 10, 1,0.1, respectively, (b) Density ratio s 
for similar values of Mo and /3. 
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Figure 2 



